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Low energy theorems are derived for the coefficients of the effective range expansion in 
s-wave nucleon-nucleon scattering valid to leading order in an expansion in which both m, r 
and 1/a (where a is the scattering length) are treated as small mass scales. Comparisons 
with phase shift data, however, reveal a pattern of gross violations of the theorems for all 
coefficients in both the 1 So and 3 Si channels. Analogous theorems are developed for the 
energy dependence e parameter which describes 3 Si~ 3 Di mixing. These theorems are also 
violated. These failures strongly suggest that the physical value of m n is too large for the 
chiral expansion to be valid in this context. Comparisons of m n with phenomenological 
scales known to arise in the two-nucleon problem support this conjecture. 


I. INTRODUCTION 


Chiral perturbation theory (yPT) has proven to be an exceptionally powerful tool for analyzing interactions 
involving pions in hadronic physics. As noted by Weinberg , however, the implementation of conventional 
yPT in low energy nuclear physics is problematic since the s-wave scattering lengths are important distance 
scales which are much larger than . This “unnatural scattering length” problem has received considerable 
attention during the past several years |i 271. The goal of this work is to implement the ideas of effective 
field theory (EFT) |2^] to a problem in which both m, and 1/a are light scales. Much of the power of EFTs 
stems from the fact that the calculation of observables is formulated in terms of an expansion parameter so 
that one has an a priori estimate of the accuracy of the calculation. However ,many of the attempts to deal 


with the “unnatural scattering length” problem, beginning with Weinberg’s JlJ, formulate the expansion at 
the level of a potential (or two-particle irreducible kernel) rather than the physical observables. In such a 
formulation, one has no immediate estimate of the accuracy of particular observables. Recently a scheme 
was introduced in which observables can be expressed in terms of a consistent power counting scheme 
In principle, such an approach realizes the full power of the EFT idea. 

The essence of this scheme is power counting in a single scale which we denote as Q: 


m-K ~ Q 1 /a~Q k~Q (1) 

where a is an s-wave scattering length, and k an external momenta. All other scales in the problem are 
assumed to be heavy and will be denoted A. The power counting implies that the external momenta must 
be small in the sense that k/A -C 1; however, one must include all orders in ka and k/m n . The inclusion 
of ka and k/m^ to all orders introduces important nonperturbative effects. This power counting scheme 
has been implemented using dimensional regularization with the PDS scheme [H Alternatively, one may 
implement it directly in configuration space using a cutoff |26| . In the case of s-wave scattering, it was shown 
that the two implementations are completely equivalent at next-to-leading order J26| and it is plausible that 
this equivalence holds generally. We refer to any scheme which systematically implements eq. (0) in the 
calculation of observables as the Q counting scheme. A number of observables have been calculated using 
Q counting at low order: N-N scattering |£l[, deuteron electro-magnetic form factors Q and deuteron 
polarizabilities (23) . All of these calculations appear to describe the experimental data reasonably well. 
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In this paper we use Q counting to derive low energy theorems for coefficients of the effective range 
expansion (ERE). The ERE is a useful way to parameterize s-wave scattering when the scattering length is 
unnaturally long. It is given by 

kcot(S) = — — + \r e k 2 + V 2 k 4 + V 3 k 6 + V 4 k 8 + ... (2) 

a 2 

As we will show below, the ry coefficients ( i > 2 ) can be calculated in a next-to-leading order calculation 
of fccot(i5). One obtains predictions which are valid up to corrections of relative order Q/A, where A is the 
scale characterizing the short-ranged physics. Moreover, these are direct predictions with no free parameters; 
in this sense these are low energy theorems. However, when we compare these predicted iy coefficients with 
ones extracted from a fit to the scattering data, we find that they fail quite badly; they predict ly’s several 
times larger than those extracted from the data. We also derive theorems, valid to corrections of relative 
order Q/A , for coefficients characterizing the momentum dependence of 3 Sq~ 3 Dq mixing at low momentum. 
These theorems also do not work well; they predict coefficients which differ from the those extracted from 
data by more than 100%. 

This poses an interesting question. Why are these theorems failing so badly when other predictions of Q 
counting seem to work rather well? One likely reason is that the successful predictions of the Q counting 
depend on different physics than what is being tested by these low energy theorems. The present predictions 
are extremely sensitive to the role played by the pion exchange in the two nucleon s-wave channel; most of 
the high quality predictions are not. It is useful to note that the Q counting EFT program of refs. p~i] . [j(i| , 
has two main components—i) dealing with the “unnaturally long scattering length” problem (Q ~ 1/a), 
and ii) “including pions” ( Q ~ to„.) p5| . It is important to test both. In Q counting 1/a and m n are both 
formally treated as being of order Q. In reality, however, m v 3> 1/a in both the singlet and triplet channels. 
Accordingly it is possible that most of the success seen to date is coming about because the expansion in 
1/a is working even if the expansion in m n is failing. 

The importance of testing the “including pions” part of the program should be very clear. Much of the 
original motivation for developing the effective field theory approach for nuclear scattering was to exploit 
chiral symmetry and develop a description of nuclear phenomenon in terms of a controlled chiral expansion. 
The advantage of including pions as explicit light degrees of freedom, rather than integrating them out as 
heavy, should be obvious. If the chiral expansion is under control, the inclusion of explicit pionic degrees of 
freedom will significantly improve the predictions at low k and will significantly increase the maximum value 
of k for which the effective treatment is useful. This suggests a simple way to test the “including pions” 
part of the Q counting program. One should calculate qualities in two EFTs—one including explicit pions 
and implementing the Q counting in eq. (Jlj), and the other in an EFT with pions integrated out and using 
Q ~ k, 1/a as the basis for the power counting. The differences between these two calculations (done at the 
same order in Q) is a measure of the effect of “including pions”. 

For many observables, the pion integrated out theory may give extremely accurate predictions. Both the 
quartet nucleon-deuteron scattering length and the deuteron charge radius [^4| were predicted using 
systematic EFT approaches with the pion integrated out. Clearly the two preceding successes reflect the 
usefulness of the Q ~ k, 1/a counting scheme (which is essentially the H counting introduced by van Kolck 
p6[), and have nothing to do with chiral counting. It is possible that most of the successes of the Q counting 
approach similarly do not test the pionic aspects. Accordingly before concluding that chiral perturbation 
theory aspect of Q counting in nuclear physics is under control, it is critical to find observables which test 
the chiral physics. 

The Vi coefficients in the effective range expansion are an ideal way to test the pion physics. The central 
reason for using the effective range expansion is to isolate the effects of an unnaturally small scattering length 
from the remaining terms in the expansion. Consider first a theory with pions integrated out. In such a 
theory, the Vi coefficients must be O(Q 0 ) since they are insensitive to a and there are no other light scales 
in the problem. In contrast, as we will explicitly demonstrate below, these coefficients are non-analytic in Q 
with Vj ~ Q- 2 ^ 2 . Thus, they are large in a theory with explicit pions and order unity in a theory with pions 
integrated out. Since there is a large difference between the two cases the predictions for the v/ s sensitively 
test the pion physics. We will also show below that when pions are explicitly included, the calculated Vi are 
not merely non-analytic in Q. For all of the Vi, each of the leading order terms in Q are non-analytic in 


2 


m w ] diverging as —> 0. Since these coefficients diverge as m n ^ 0 it is clear that they are dominated by 

pionic effects. 

This paper is organized as follows: In the following section we derive low energy theorems for the coefficients 
in the effective range expansion; analogous theorems for 3 S±- 3 Di mixing are derived next. After that we 
compare the predictions with coefficients extracted from the Nijmegen partial wave analysis of the 
scattering data. Finally we discuss why these failures might have been anticipated in light of the known 
scales in the nucleon-nucleon problem, and how the failures of the low energy theorems can be reconciled 
with the apparently successful predictions for the phase shifts reported in ref. [^lj . 


II. LOW ENERGY THEOREMS FOR THE EFFECTIVE RANGE EXPANSION 


In developing physical intuition for our low energy theorems, we find it somewhat more useful to use the 
cutoff formulation discussed in ref. |26| . The essential physical idea in this approach is to implement the 
separation of long distance physics from short distance physics directly in configuration space. A radius, R, 
is introduced as a matching point between long and short distance effects; renormalization group invariance 
requires that physical quantities must be independent of R. It is important, however, that R be chosen large 
enough so that essentially all of the effects of the short distance physics is contained within R. The potential 
is divided into the sum of two pieces, a short distance part which vanishes for r > R and long distance part 
which vanishes for r < R. At R, the information about short distance effects is entirely contained in the 
energy dependence of the logarithmic derivative (with respect to position) of the wave function at R. Thus, 
provided we can parameterize this information systematically, we can formulate the problem in a way which 
is insensitive to the details of the short distance part of the potential. For r > R the Schrodinger equation 
is solved subject to the boundary conditions at R. For s wave scattering, the wave function at R may be 
parameterize as Asin(fcr + <5o); the energy dependence of the logarithmic derivative is independent of A and 
can be expressed in terms of an expansion similar to an effective range expansion: 

A: cot(<S 0 ) = — l/o s hort + 1/2 r° k 2 + v$ k 4 + v° k 6 + k 8 + ... (3) 


Power counting in Q for s wave scattering can be implemented straightforwardly. All of the coefficients in the 
proceeding expansion are assumed to be order Q° except the first term (—l/a s hort ) which will be taken to 
be order Q 1 to reflect the unnaturally large scale of the scattering length. Power counting for the long range 
part of the potential simply follows Wienberg’s analysis (lj, with the previso that the potentials are only used 
for r > R. At order Q 2 in fccot((S), only the simple one pion exchange contribution to the Vi ong contributes. 
The power counting also justifies an iterative solution of the Sclirdinger equation for r > R along the lines 
of a conventional Born series. It differs from the usual Born series in that the boundary conditions at R are 
implemented. Finally, Q counting is used in expanding out the final expression for fccot(<5). 

Carrying out this program gives the following expression for fccot(<5) at order Q 2 for the 1 Sq channel 
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We use the convention in which f n = 93 MeV. Apart from well-known parameters from pionic physics, there 
are three parameters— cio , d, and r°. Where and l/a s hort from eq. (|./j) is rewritten as 1/ao + dm 2 , with 
1/ao ~ Q and dm 2 ~ Q 2 . These parameters fix the energy dependence of the wave function at the matching 
scale i?; renormalization group invariance requires d to depend on R logarithmically. The parameter ao 
corresponds to the scattering length at lowest order in the Q expansion. It is related to the observed 
scattering length by 
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The PDS scheme is superficially quite different in mathematical detail but ultimately describes the same 
physics based on Q counting. In terms of the parameters in the PDS formulation, one obtains the same 
expression [^ 6 | provided one relates the Co, D 2 and C 2 coefficients in terms of do, d and r° according to 
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There is a subtle issue associated with the equivalence given above concerning the behavior as the chiral 
limit is approached. This issue is discussed in the appendix. 

The parameter do plays a critical role in eq. ([|). It not only plays a dominant role in fixing the value of 
k cot (5) as k —> 0 , it also appears in the tan -1 term and one of the In terms which have nontrivial dependence 
on k/m n . Unfortunately, we do not have a direct experimental way to fix do- On the other hand, from eq. (^) 
we see that 1/ao = l/a(l + 0{Q/ A)). Accordingly, if we replace 1/do with 1/a in the terms of order Q 2 in 
eq. (|]), any error made is order Q 3 —which is beyond the order to which we work. Thus, to order Q 2 , 
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The form of eq. (Jt[) also holds for the triplet channel—except for changes in the value of the parameters. 
The distinction between the singlet and triplet channel is only due to the tensor force. The effect of the 
tensor force on the s-wave requires at least two iterations of the pion exchange and contributes to k cot(5) 
only at order Q 3 and beyond pT| . 

The predicted coefficients in the effective range expansion are easily obtained. One begins with the 
expression for k cot((5) in eq. (i 7 j) and simply makes a Taylor expansion with respect to k 2 . From the definition 
of the effective range expansion in eq. (Ejj) the 77 coefficient is simply the coefficient in the Taylor expansion 
multiplying k 2i (for i > 2) We find: 
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The expressions for the 77 coefficients in eq. (||) are low energy theorems. As is clear from our derivation, 
they are valid provided Q/A is sufficiently small. As stated in the Introduction, one sees that 
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Moreover, it is clear that each term in the equations for all of the terms in the expressions for the 77 go as 
1 /to” with n> 2. Thus the Vj all diverge as the chiral limit is approached. 
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III. LOW ENERGY THEOREMS FOR 3 S, AND THE 3 D 1 MIXING 


We can also derive low energy theorems which test the pionic contributions in the mixing between 3 Si 
and the 3 D i. Conventionally the S matrix in this mixed channel is parameterized in the form 

/ cos(2e) sin(2e) , 

* “ V sin(2e) cos(2e) e i{2S °> ) ^ U ' ) 

The leading order contribution |pl|| to the e parameter in Q counting is 0(Q 4 ); when pions are integrated 
out, however, the leading contribution is 0(Q 2 ). Thus the e parameter provides a sensitive test of the pion 
physics, in the sense discussed in the Introduction. 

Rather than consider the e parameter globally as a function of fc, it is illuminating to expand the e as a 
function of k: 


zip) = 3i k 3 + 32 k 5 + g 3 k 7 + ... (11) 

In a manner analogous to the derivation of the theorems for Vi we derive theorems for the 3,; coefficients: 
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As with the Vj, each term in the leading order expression for gj diverges in the chiral limit of m, r —> 0. 


IV. COMPARISON WITH SCATTERING DATA 

The Vj and gj coefficients are observables which can be extracted from the nucleon-nucleon scattering 
data. We have extracted these coefficients from the scattering data as parameterized in the Nijmegen fl29|] 
partial wave analysis. The extracted coefficients for the Vj in the triplet channel are taken from ref. [ [l0[ . 
For the Vj in the singlet channel we have done a least squares fit directly from the Nijmegen phase shifts at 
very low energies. Our fits in the singlet channel agree quite well with fits of the coefficients calculated from 
potential models fitted to the phase shift data ,0- We have also extracted the g t coefficients using a least 
squares fit to the low energy data. In Table (IJ) we compare the extracted coefficients with those predicted 
by the low energy theorems. 

It is apparent from Table (|) that the low energy theorems for the u, coefficients fail quite badly in both 
the singlet and triplet channels. In all cases they predict coefficients which are several times the extracted 
ones, typically by a factor of ~ 5. The V 2 in the triplet channel has a more spectacular failure—overstating 
the extracted value by more than a factor of 20. The low energy theorems for the gi coefficients also do 
rather poorly. The predicted gi coefficients are all more than 100% greater than those extracted from the 
data, essentially meaning no prediction. The coefficient 32 is ~ 4 times that extracted from the data while 
predicted 33 is ~ 8 the extracted one. 
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V. DISCUSSION 


The clear failure of the chiral expansion to describe the Vi coefficients for the s-wave may seem at first 
glance quite surprising. After all, Kaplan, Savage and Wise described s-wave scattering using the PDS 
approach and, using a global fit for the parameters, appear to successfully describe the phase shifts up 
to k ~ 300 MeV. How, then, can the Vi coefficients in the effective range expansion describing the same 
scattering data all be badly wrong? The answer is, in fact, quite simple. The principal point is simply that 
the general shape of the phase shift curve—a very rapid rise at low k followed by a slow decrease is implied 
directly from the effective range expansion with large a; any theory with free parameters consistent with the 
effective range expansion will be able to reproduce the crude shape. Thus the ability to crudely reproduce 
the shape of the data is not a stringent test. Moreover, the global fit used in ref. J|l]] masks possible subtle 
problems with the detailed shape of the fit to the phase shifts. However, there is an important issue, related 
to the quality of the fit of the s-wave scattering in ref. |2lj] which goes to the heart of the issue relating to 
the role of pions in the effective field theory program. 

As we have argued in the Introduction, the easiest way to test the role of pions in the EFT is to compare 
a calculation in an EFT with explicit pions to a calculation at the same order in an EFT with the pions 
integrated out. If the chiral expansion is working well we should find that the inclusion of explicit pions 
both substantially improves the quality of agreement of the low k observables with the data and increase 
the range in k for which the EFT is useful. Here, since we are working at order Q 2 in kcot(S ), the theory 
with pions integrated out is simply the effective range expansion up to the effective range term. 

In fig. 0 we compare the difference of fccot(<5) from the Nijmegen j2l| partial wave analysis in the 1 S'o 
channel from three order Q 2 EFT predictions: a) a pions-integrated-out theory with parameters fit to 
the observed a and r e (this is precisely the first two terms in the effective range expansion); b) a pions- 
explicitly-included theory with the parameters in eq. (iTj) fit to the observed a and r e using eq. (||); c) a 
pions-explicitly-included theory from the global fit of Kaplan, Savage and Wise ©• 

We have plotted the difference from the data to focus attention on the accuracy of the predictions (in 
the spirit of refs. [ ^o|j27| ). Plotting the data in this manner also removes the visual effect of reproducing 
the rise and fall of the phase shifts. The striking thing about fig. |I] is that the inclusion of pions does 
not improve the quality of the fit at low k compared to the pion-less theory. Indeed, at moderately low 
energies, the inclusion of pions markedly worsens the agreement of theory with experiment compared to the 
pion-integrated-out theory. This strongly suggests that the chiral expansion is not under control. Having 
seen this, it is not so surprising that the global fit is so poor at low energies; it was designed to compromise 
the low energy behavior in order to do better at larger Q. In testing low energy observables, however, it is 
more appropriate to use coefficients fit to the low energy behavior. Indeed, there is an ambiguity in using 
the global fit parameters. In the plot we have evaluated k times the cotangent of the phase shifts from the 
expression in ref. [;2lj . Alternatively, we could have used the global fit parameters directly in the expression 
for k cot <5 in eqs. and ([|). While these two forms agree to order Q 2 , they differ at higher order. As the 
fit was done using 6 rather than kcot(S) we felt it more reasonable to plot it this way. Had we used the 
same parameters in the fccot(d) expression the fit would have markedly worsened. 

The Vi coefficients directly measure the improvement of the theory with explicit pions over the theory 
with pions integrated out. The pion-less theory at this order is simply the effective range expansion up to 
the r e term; the v % coefficients parameterize the way in which the phase shifts deviate from this expression. 
Thus the failure of the effective field theory treatment to get the coefficients is also evident in fig. [I] and 
suggests a failure of the “including pions” part of the Q counting program—at least so far as the s-wave 
scattering is concerned. 

It should also be remarked that the failure in the prediction of the e parameter evident in these gt 
coefficients is much worse than suggested by the plot of e against k in ref. |nj. In ref. , at low fc, the 
predicted e appears to be ~ 40% above the data at low energy (although the predicted shape differs greatly 
from the observed one at higher energies). How can this be reconciled with a g\ parameter predicted in our 
low energy theorem which is 2.3 times the observed one? The origin of this discrepancy is related to the fit 
used. Using the global fit parameters of ref. (2l) and directly calculating g\ yields a coefficient which is only 
~ 40% greater than that extracted from the data. In contrast, using the physical scattering length, as we 
do in eq. we obtain one which is ~ 130 % larger. As both the calculation using the physical a and the 
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global fit are at the same order in Q , one must view the better agreement of the calculation using the global 
fit parameters as being fortuitous. Moreover, the large disagreement between calculation based on the global 
fit and the one based on the a is itself an indication that the expansion has broken down. 

It should be obvious from the pattern of failure for both of these low energy theorems that the A, the 
scale of the short distance physics, is not much larger than m n . The derivation in ref. p6| provides some 
insight as to why this might be expected which may not be as immediately apparent in the PDS formalism. 
In the derivation using the cutoff formalism, the expansion in m, r /A came in two distinct ways. The first 
was in justifying the perturbative inclusion of the pion exchange interaction in a manner similar to the Born 
series but with boundary conditions at some matching radius R to reflect the short distance physics. The 
second was in the expansion of m^R in the various integrals which arise. For this expansion to be valid 
we must take m^R « 1 . A central issue is whether this is true in practice. Of course, the value of R is 
somewhat arbitrary having been introduced by us. We are not free, however, to make R arbitrarily small; by 
construction, R must be large enough so that the effects of the short distance potential on the wave function 
are contained within R with an accuracy comparable to the order at which we are working. 

In fact, comparisons with potential models suggest that m^R ~ 1, i.e. the range of the short distance 
potential is comparable to 1/m n . One rather compelling way to see this is found in sect. IV of Scaldeferri 
et al Q. There, using a generalization of conventional effective range theory, it is shown that for any energy 
independent potential model consisting of a one pion exchange plus a short distance potential (including a 
nonlocal potential), there is a rigorous lower bound on R. In particular, it is shown that in order to fit the 
scattering length and effective range, the short range potential must make a non-negligible contribution to 
the wave function out to at least 1.1 fm— i.e. R > 1.1 fm. This corresponds to m^R > .78. Moreover, it 
is clear from the derivation that this bound on R cannot be saturated, since the bound is saturated only if 
the wave function is strictly zero for r < R; the actual value of R beyond which the short distance potential 
ceases to make significant contributions is accordingly expected to be significantly larger suggesting a value of 
m^R ~ 1. For example, if the short distance potential for the 1 S'o channel is taken to be a square well whose 
depth and range is fit to the scattering length and effective range, then R = 2.3 fm ||], which corresponds to 
m^R « 1.6. An obvious conclusion is that the scale of the short distance contribution to the nucleon-nucleon 
interaction is not well separated from m w making the chiral expansion highly problematic. 

Of course, it might be argued that the derivation in ref. ||] includes only one pion exchange and short 
distance physics. It does not include higher order long distance effects such as two pion exchange which 
could alter the result. However, if one is in regime where Q counting is valid, then two pion exchange and 
similar mechanisms only contribute to kcot{5) at order Q 3 or higher and may be neglected, implying that 
the bound in ref. || is valid up to corrections of order Q/A. 

The justification of the perturbative, Born-like expansion for the inclusion of the pion exchange interaction 
also may be problematic. In ref. |21[] the quantity, 16nf 2 /(g\M), was taken to play the role of A and 
the dimensionless parameter, r] = (m^g^M) / (I6nf 2 ), was identified as playing the role of the expansion 
parameter. In terms of the derivation of ref. (g 6 ), it is clear that 77 = fi(k 2 = 0). The potential difficulty for 
X PT was identified in ref. pl[ 77 is rather large; numerically 77 « .47. If the relevant parameter is in fact 
77 , one might expect a slowly converging theory. On the other hand, if the relevant parameter is more like 
2 77 , things are quite out of control. Moreover, this expansion parameter arises in the context of the central 
force. In the triplet channel there is also a tensor force. While formally in Q counting the effect of the tensor 
force is higher order, experience in the nuclear physics problem is that it is quite important owing to large 
numerical factors. Given both the problems with the m^R expansion and the justification of a perturbative 
treatment for the pion, the failures of the low energy theorems are not so surprising. 

The possibility of problems with the perturbative treatment of the pion exchange interaction implicit in the 
PDS scheme has been raised previously by Gegelia |l9| and Steele and Furnstahl |27| using rather different 
arguments. Here we have shown explicitly in the case of the 1 4 coefficients in s-wave scattering and the gi 
coefficients in the 3 Si- 3 D\ mixing that the “including pions” aspect of the systematic Q counting approach 
in nuclear physics fails badly. We have also given general arguments that there exist “short distance” scales 
in the nucleon-nucleon problem which are ~ ra r suggesting that the “including pion” parts of Q counting 
may be expected to have problems. It remains to be seen just how widespread these problems will turn out 
to be. Clearly the perturbative pions work well in high partial waves, but the question remains whether 
observables sensitive to both the pion physics and to s-wave nucleon-nucleon interactions are generally well 
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described in Q counting. In order to determine this, it is necessary to study a number of observables which 
are particularly sensitive to pion range physics. 

The authors thank Silas Beane and Daniel Phillips for interesting discussions. We are indebted to Vincent 
Stoks for providing us with numerical values of the effective range expansion coefficients for potential model 
calculations of the 1 S'o channel. We also thank Mike Birse for pointing out to us the issue discussed in the 
appendix. TDC gratefully acknowledges the support of the U.S. D.O.E. 


APPENDIX A: THE CHIRAL LIMIT AND THE EQUIVALENCE OF THE CUTOFF AND PDS 

APPROACHES 


In ref. (26) it was shown that the PDS and cutoff expressions for k cot S are equivalent provided the 
coefficients in the two approaches are matched as given by eq. (§• However, as pointed out by Birse |3^] 
there is an apparent difficulty with this equivalence as the chiral limit is approached. Note that Q counting 
formally holds if both 1/a and m^ are small, however the combination m^a is unconstrained and can take 
any value. Thus, a scheme which implements Q consistently should be able to correctly describe the formal 
limit TOtt —s> 0, 1/a —> 0, m^a —> 0. This formal limit, is of course nothing but the chiral limit in a regime 
where Q counting holds. We note, at the outset, that this limit is clearly not of immediate practical concern 
since in nature m^a 3> 1. Never-the-less, if the Q counting scheme is viable, nothing in the formalism should 
prevent an approach to the chiral limit. The potential problem is quite clear in the last equation of eqs. ®): 
the left-hand side goes to zero in the chiral limit while the right-hand side does not, suggesting that the 
equivalence cannot hold for all values of m n . 

The difficulty stems from a minor inconsistency with the PDS formulation which can be cured easily 
and which has no observable consequences for NN scattering. Note that as the chiral limit of the last 
equation of eqs. ()(]) is taken, the right-hand side is not merely nonzero, it is /i dependent. This is a direct 
consequence of the fact that in the formulation of PDS in ref. (2l|], the scattering length in the chiral limit 
is not renormalization group invariant. The PDS expression for kcot(S) at order Q 2 is J2(| 
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Taking the limit k —> 0 picks up — 1/a and subsequently taking the chiral limit yields: 
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which depends on g and thus is manifestly not renormalization group invariant. Note, this problem is not 
immediately apparent in ref. |2l) since the expression for the scattering length is only given after /i is set 
equal to m n . 

The reason for the g, dependence of eq. (A2) is clear. No counter term has been given in ref. [^l| which can 
absorb it. On the other hand, this can be dealt with easily. The simplest way to rewrite the coefficient Co 
as c£ on P ertub _|_ Qpeiturb w j iere ^nonpertub j teratec j to or ders and contributes to k cot(<5) at order Q and 

reproduces the physics 1/ao- In contrast (jP cltulh serves as a counter term absorbing the g dependence, is 
not iterated to all orders and only contributes at order Q 2 . Note, however that in calculating cross-sections 
^perturb an j Dm 2 enter in exactly the same way. Thus, only the combinaton ( 7 P erturb _(_ Drn 2 ^ affects the 
scattering amplitude. To the extent that D was fit directly from scattering data one can include (7P elturb 
by making the substitution D^rn^ — > Dm J 4- QP ert " rb w ith no change to the predicted scattering. Thus, 














the only effect of including £;P eltulb as a separate counter term is on the behavior as the chiral limit is 
approached. 

After the work in this appendix was finished, we became aware of the work of Mehen and Stewart 
in which substantially the same conclusions are reached, namely that /i independence requires splitting 
a perturbative contribtuion from Co and that only the contribution from the combination of the D plus 
the perturbative part of Co contributes in nucleon-nucleon scattering. We note that the work of Mehen 
and Stewart is far more general than this and proposes a more transparent way to implement dimensional 
regularization than PDS. 
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5 ( 1 5'o channel) 

V 2 (fm 3 ) 

V 3 (fm 5 ) 

V 4 (fm 7 ) 

low energy theorem 

-3.3 

17.8 

-108. 

partial wave analysis 

-.48 

3.8 

-17. 


5 ( 3 Si channel) 

V 2 (fm 3 ) 

V 3 (fm 5 ) 

V 4 (fm 7 ) 

low energy theorem 

-.95 

4.6 

-25. 

partial wave analysis 

.04 

.67 

-4.0 


e ( 3 Sj- 3 Di mixing) 

3i (fm 3 ) 

32 (fm 5 ) 

33 (fm 7 ) 

low energy theorem 

3.9 

-86. 

1.8 10 3 

partial wave analysis 

1.7 

-26. 

2.2 10 2 


TABLE I. A comparison of the predicted effective range expansion coefficients, Vi, for the 1 So and 3 Si channels 
and the predicted coefficients in the expansion of e with coefficients extracted from the Nijmegen partial wave 
analysis. 
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FIG. 1. k cot($EFT) — k cot(^NPWA) versus k where the subscript NPWA indicates the Nijmegen partial wave 
analysis and the subscript EFT indicates the effective field theory. 
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